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ABSTRACT 


This project is an effort towards modelling and 
developing an electrical simulation of various activities 
associated with the nervous and cardiovascular systems. 

The neural cell model has been simulated by a circuit 
incorporating its various analog and digital properties , while 
disregarding the shape of action potential waveform. On the 
other hand, a simulation of the membrane patch with its 
passive and active properties and an approximation of action 
potential waveform has been developed. 

Using the neural cell model as a basic module, 
peripheral nervous system and a spinal reflex have been 
simulated. Higher activities of the central nervous system 
have been simulated on a microprocessor-based system, 
incorporating six sensory and four effector channels. By 
changing a part of software, different kinds of conditioned 
reflexes* learning processes and vegetative activities can 
be simulated . 

A simplified simulation of blood circulation system 
alongwith the electrical and pumping activities of the heart 
has been done. In addition to, automatic mechanism of the 
heart itself, vegetative control by the central nervous syster 
also has been simulated. Effects of bleeding, transfusion, 
physical activity and disorders in the conduction of excitatic 
can be studied on it. Through software, effect of external 
stimuli on cardiac activities can also be incorporated. 
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CHAPTER I 


• INTRODUCTION 

The problem of the relation between animate and 
inanimate matter has excited the interest of mankind since 
ancient times. Investigations in this direction have proceeded 
along two lines. The first of these aims at the artificial 
creation of living matter. Despite the efforts in this 
direction throughout the history of mankind, the aim has not 
yet been achieved. The efforts in the second line have 
attempted to simulate various functions of living organisms by 
means of inanimate systems. Studies and investigations 
associated with developing the engineering models simulating 
the structure or behaviour or 'both of living organisms 
fulfil several objectives. They help in learning and 
understanding biological problems. The understanding developed 
in the process of developing such analogs may lead one to 
design various kinds of engineering devices utilising to a 
greater or lesser extent, the principles of the structure of 
living organisms. 

1.1 Modelling Techniques 

In developing analogs, both the mathematical and 
physical modelling methods have been used. Mathematical 
models have great utility in limited domains. They are 
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invaluable in cases where the number of variables is reasonably 
limited and nonlinearities do not present severe analytic 
difficulties. They are of vital importance, when statistical 
or probabilistic aspects are decisive such as in the membrane 
biophysics or the neural network formulations. 

Physical models are relatively more tangible than their 
mathematical counterparts, depending on how they are made, 
they may be placed in one of following catagories 

(i) Ionic models 

(ii) Pneumatic and electro-mechanical model 

(iii) Electronic models 

(iv) Computer simulated models 

1.1.1 Ionic models 

These are based on the effects associated with electro- 
chemical reactions. The most important application of this 
approach has been in modelling the processes in a nerve cell. 
However, this approach does not find an extensive use, because 
the measurements may turn out to be as difficult as in the 
real situation and further the physical mechanisms may turn 
out to be as obscure as the real phenomena. 

1.1.2 Pneumatic and electromechanical models 

The most extensive use of such models has been in 
simulating blood circulation system and physical movements. 
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They have been also used for simulating some features of the 
cell such as adaptational aspects. However because of .their 
cumbersome nature, they are less common. 

1.1.3 Electronic analogs 

These can simulate continuous variable nonlinear 
operations economically and quite often accurately. They 
permit a rapid and effective kind of observer-mod'el interaction 
as the different time-dependent phenomena can be directly 
observed while stimuli and model parameters are changed. 

However, electronic models are not good solution for simulating 
learning and other brain processes envolving complex neural 
networks, because the observation and manipulation of parameters 
become very difficult. 

1.1.4 Computer simulated models 

My computer simulated model is essentially an extension 
of some mathematical model since a set of mathematical 
equations describing the given phenomena is generally the 
starting point for this kind of modelling. In this approach, 
either the analog or the digital or both kind of computers are 
used, inalog computers permit the direct observation of 
phenomena, but usually the models based on them tend to be 
cumbersome as a large no. of modules are needed. The digital 
processing intrinsically does not permit economical 
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representation of continuous variable nonlinear interactions. 
However, the special problems that arise in large network 
simulation are more readily handled by digital computers than 
by other techniques. Due to the growing speed and storage 
capabiliti es , they may ultimately provide the most satisfactory 
means for modelling complex neural systems. 

1.2 Scope of the present work 

Most challenging problem in the area of simulation of 
structure and behaviour of living organisms has been the one 
associated with man himself. As the human body involves a 
tremendously complex set of processes, in any attempt towards 
simulating one has to limit the scope of the work. In the 
present work, an effort has been made in the direction of 
simulating some very simple processes related with the nervous 
system activities. As an example of various processes in the 
body that are controlled by the nervous system, as well as 
have their own automation mechanism, the cardiovascular system 
has been chosen. 

In the second chapter, the properties and features of 
nerve cell - the basis of entire nervous aysteci are presented. In 
third chapter, models simulating nerve cell have been proposed. 
In the next two chapter, the organisation and function of the 
nervous and cardiovascular systems are d iscussed. After 
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presenting a picture of these systems and developing nerve 
cell models, the simulations developed for the cardiovascular 
and nervous systems are described in the sixth and seventh 
chapters respectively. A discription of the overall system 
is presented in the chapter eight. ThtS concluding chapter 
suggests the experiments and investigations that can be carried 
out using the system developed in the project. 



CHAPTER II 


NERVE AND MUSCLE CELLS 

A cell is a basic unit of life capable of independent 
manifestations of life. In living organisms, it is possible 
to find cells widely differing from one another. This is the 
result of the process of evolution necessiating specialisation 
of individual cells. Aggregates of cells that perform a 
common physiological function are called tissues. These can 
be divided into four classes : 

(a) Various kinds of epithelia - They form the basis of 
glandular control. 

(b) Muscular tissues - They consist of elongated cells 
that are capable of contraction when excited by stimuli and 
thus convert chemical energy directly into mechanical energy. 
The action of muscular tissue is mostly controlled by the 
nervous system. 

(c) Connective tissues - They bind together the other 
kind of tissues. 

(d) Nervous tissues - They form the basis of nervous 
control. 

Despite the considerable differentiation of all cells.', s 
there exist certain properties common to all. However, the 
nerve and muscle cells have special property that a temporary 
change in trans-membrane potential (called action potential) 
is developed upon application of excitation. 


7 


2.1 Nerve Cells 

A nerve cell (also called neuron) is the fundamental 
element of the nervous system. Individual nerve cells may 
differ from one another in anatomical configurations and in 
functional properties. However, all neurons, regardless of 
individual differences in shape, size and location have four 
functional components : 

(a) one or more input elements 

(h) an integrative element 

(c) an active transmission line 

(d) one or more output elements 

The integrative element is soma or the cell body. There 
are nerve processes that grow out of the cell body proper and 
may have a large variety of shapes. The cell body is similar 
in structure to the cells of' other tissues and its diameter 
varies from several microns to several tens of microns. The 
processes are small in diameter while their length may run 
to anything from a fraction of a millimetre to the order of a 
metre. The inputs to a neuron may occur at several points 
on its surface. The majority of the inputs, however, enter 
through short, highly branched processes, called dendrites, 
which at their terminal ends interface with extensions from > 
other neurons or sensory cells. A long process of 
approximately uniform diameter is called on axon, which acts 
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as the active transmission line. The end of the axon is 
branched to form terminal arborization. 

The integrative process either takes place in the 
dendrites or in the soma. If the sum of input excitations 
exceeds a threshold level, the cell fires, generating a signal 
•which is actively transmitted down the axon length to the 
terminal regions. The transmission properties of the axon 
are bidirectional. However, the nerve cells are connected to 
one another by means of special nerve endings, called synapses. 
They Incorporate an electrochemical mechanism for information 
flow from one cell to other cell in one direction only. 

The nerve cell like every cell in a living organism 
contains a fluid, called the intracellular fluid (or the cell 
plasm), which is separated from the interstitial fluid 
surrounding and bathing the cell by an extremely thin 
(50-100 A.U. ) membrane. Both of the fluids are ionic 
solutions, the predominant ions in the interstitial fluid 
being Ua + and Cl”, while the intracellular fluid consists 
essentially of K + ions and an organic ion generally referred 
to as A”. The electrochemical processes which give rise to 
most of the properties of the nerve cells are controlled 
principally by the permeability of the membrane to these ions. 
In the so called myelinated axons, the membrane is covered 
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"by a much thicker (2-3 /il) lipo-protein sheath, thereby 
rendering the membrane impermeable to these ions. The myelin 
sheath is broken at regular intervals (0.4-2 mm) for about , 
thus allowing the migration of ions across the membrane only 
at these regions, called the nodes of Ranvier. 

2.2 Action Potential 

The membranes of nerves and some other cells posses 

the highly distinctive property of being excitable. A potential 

difference exists across the membrane. This is called 

membrane potential. In the steady state, the cell remains 

polarised with the intracellular fluid at a negative potential 

of 60-100 mV with respect to the interstitial fluid. This 

potential is called the resting potential. Complex, but 

shortlived, electrochemical processes propogated on the 

membrane surface are initiated under the influence of 

particular stimuli. The electrical response during these 

[l 2 1 

processes is called action potential * . 

An environmental change brings about a transient 
depolarisation, usually by increasing the permeability to 
Ha + ions. , Influx of Na + depolarises the membrane, thereby in 
turn increasing the permeability to Ha + which leads to further 
depolarisation. If the original stimulus is capable of 
depolarising the membrane beyond a threshold value, the process 
is regenerative. Hot only the membrane gets fully depolarised 
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but even its sign changes (going about 30-70 mv). This rise 
of .impulse constitutes the first phase of action potential. 
The stimulus serves merely to trigger the cell into activity 
and an additional increase in stimulus intensity beyond 
threshold has no effect on the amplitude of action potential. 
For this reason, it is said to have all-or-none behavior. 
Towards the end of first phase, the depolarisation peak is 
followed by a fall in sodium permeability and a rise' in 
potassium permeability. This starts repolarisation process. 

Hh 

The accompanied voltage change results in a decrease in Fa 
permeability and after a delay in K + permeability also. When 
repolarisation is complete, K + permeability is still above 
resting value and consequently, the membrane hyperpolarises, 
going even more negative than the transmembrane potential* 

After the end of impulse proper, in the third phase, 
slower processes related to the membrane returning to the 
initial state occur. Fig. 2.1 depicts a typical action 
potential waveform and associated processes. 

In a short interval of time after the onset of action 
potential, the membrane is non-excitable, This interval is 
called the absolute refractory period. It specifies t-he 
maximum impulse repetition rate, usually of the order of 1000 
impulses per second. 
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After an absolute refractory period, it is possible 

but difficult to generate an action potential.. This is called 

the relative refractory period. However, in time the threshold 

diminishes to the normal value. A typical simplified picture 

~?3l 

of cell sensitivity after an excitation is given in Pig. 2.2 L J 

The plot of changes of excitability as a function of 
time implies that the higher the excitation potential (Uex) 
applied to the cell, the shorter the time between individual 
action potentials. It thus follows that the nerve cell is 
an element capable of frequency modulation. Shorter interval 
between pulses and associated higher impulse rates correspond 
to more intense stimulation. 

If the excitation stimulus be applied in the form of | 
a pulse, the strength required to initiate an impulse depends 
on the duration of the excitation pulse, the relationship is 
given as i 

1st = I rh (1 - e~ T/Z T 1 (2.1) 

where 1st = stimulus strength 

T = stimulus pulse width 
H = a parameter characteristic of the 
membrane passive properties 
I h = stimulus strength to initiate an 

impulse when . This is known as 

rheobasic stimulus. 
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2.3 Nerve Conduction 

For subthreshold stimulation, the propagation along the 
axon length is analogous to that in a passive electric cable 
and is referred to as electrotonic conduction. The velocity 
and attenuation are decided by its passive properties such as 
resistances of the intracellular and interstitial fluids, 
dielectric constant and leakage conductance of membrane etc. 

The important mode of propagation as regards the 
transmission of information is the suprathreshold propagation. 
As the stimulus exceeds the threshold, an action potential is 
developed. The process involves a short section of t,he 
membrane. However, the sudden change in the transmembrane 
potential during the first phase leads to the flow of surface 
currents and thus results in a rise in transmembrane potential 
at neighbouring points of the cell membrane. When this change 
exceeds threshold, action potentials are generated. As the 
action potential propagates along the active membrane, it 
leaves the points through which it has passed 3n refractory 
period and they will not respond to the surface currents due 
to neighbouring regions. Propagation of action potential thus 
takes place only in the direction of the incident wave. 

In unmyelinated fibres, the stimuli propagate relatively 
slowly (about 1 metre per second or so). Myelination greatly 
increases the conduction speed and thus decreases the reaction 
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time in the nervous system. The action potential is generated 
only at the nodes of Ranvier, the propagation being 
electrotonic in the myelin -sheathed segments. As this region 
has far less capacitance and leakage conductance compared to 
an unmyelinated fibre of same length and diameter, the speed 
of conduction is much higher. Owing to the attenuation 
associated with electrotonic propagation, the amplitude of the 
impulse is considerably reduced when it reaches the next node, 
but is still large enough to excite it. Thus the 
depolarization jumps from one node to the other and this mode 
of conduction is referred to as saltatory conduction or step 
conduction. 

2.4 Synaptic Conduction 

The conduction between mutually interactive nerve cells 
is through a synapse. A single axon may form synaptic 
connections to many nerve cells. The axonal side of the 
synapse is referred to as presynaptic region, the corresponding 
portion of the succeeding cell is termed as post-synaptic. 
The separation between pre-and post-synaptic membranes, called 
synaptic gap is of the order of 200 A.U. 

Under the influence of impulses arriving at the 
presynaptic knob, a special activiting substance called 
mediator is released from the vesicles at the membrane. This 
secretion induces a rise in sodium permeability of post synaptic 
menbrane resulting in an increase in transmembrane potential. 
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This is referred to as excitatory post synaptic potential 
(EPSP). At some functions* an inhibitory mediator is released, 
causing hyperpolari&ation of the membrane. This is called 
inhibitory postsynaptic potential (IPSP). 

The synaptic space also contains a substance causing 
decomposition of the mediator. Due to its action membrane 
potential returns to its original value. However if the 
postsynaptic potential exceeds the threshold, regenerative 
depolarisation starts and action potential is generated. The 
process of excitation or inhibition of the qynapse delays the 
transmission of signals from one neuron to the other. This 
delay, called synaptic delay, is of the order of 0.5 msec. 

The axonal delay, which depends on the speed of propagation 
and the length of the axon is in the range of 0.1 - 20 msec. 

A nerve cell is induced by stimuli from many synapses 
on the cell body and dendrites alike. This is called spatial 
summation. Further, since the mediator secreted by the 
synapse decomposes with a certain time constant (2-4 msec), 
impulses which arrived at the synapse earlier also have some 
contribution in the formation of net excitatory potential. 

This phenomenon is called temporal summation. 

2.5 Muscle Cells 

The electrical properties of voluntary or skeletal 
muscle cells are comparable to those of nonmyelinated axons. 
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Differences between muscles and nerves lie mainly in the 
special anatomical features of muscle cell and in the role 
played by Ca ++ ions in muscle excitation ^ "l An action 
potential initiated at one end of a muscle fibre spreads by 
virtue of local current effects to the other end. The action 
potential is accompanied by a contraction. 

A skeletal muscle is built-up from a set of individual 
10-100 ;u-thick fibres which are functionally separable. 
Increasing the stimulus strength to a muscle produces an 
increasing response until a maximum level is reached. The 
increasing response results from an increase in the number of 
activated fibres. When a stimulus initials a propagated action 
potential on a fibre, resulting contraction throughout the 
fibre is maximum. When a stimulus does not initiate 
propagation, the contraction is confined to activated region 
alone. 

Efferent (or motor) nerve fibres convey electrical 
activity to the appropriate muscle leading to its contraction.' 
Each axon normally activates from 3 to 150 muscle fibres. The 
interface between a motor nerve ending and a muscle fibre is 
known as the neuromuscular junction. The mechanism involved 
is similar to the one in synapse. The neuromuscular delay 
ranges from 0.5 to 1.5 msec. 
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2.6 Receptors 

Physical stimuli like light, heat, pressure, odour 
etc are converted into neural signals by a set of specialised 
cells known as receptors. Stimulation of a receptor cell 
produces a decrease in membrane polarisation, called the 
generator or receptor potential. Its value is related to 
the intensity of the stimulus and is graded as well as 
stationary. Each receptor is primarily responsive to a 
particular form of energy % thus there are optoreceptors, 
presso-receptors, chemoreceptors etc. Some receptors 
respond only to unidirectional change while some are 
primarily responsive to the rate of change in the stimulus 
intensity. 

If a steady state stimulus is applied, the resulting 
generator potential will ordinarily diminish with time. 

This decrease in potential is called adaptation. For some 
kind of receptors a decrease from the initial value to a 
lower steady state value is observed.^ , ^ ,r ^ 
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NERVE CEIL SIMULATION 

A nerve cell is an analog digital element, and is thus 
distinguished from the fundamental components used in data 
processing and automatic control systems. The following is a 
list of well established properties of the nerve cell^ ^ ~ 

i) A cell may have a number of excitatory as well as 
inhibitory input channels and only one output channel, 
viz - axon. 

ii) The input channels (through synapses or dendrites) 
incorporate spatial and temporal summation. 

iii) The output of a nerve cell follows the " all or none " 
principle. 

iv) in impulse or a series of impulses is generated when 
the net excitation exceeds a certain threshold. 

v) Post synaptic changes in excitability (absolute and 
relative refractory periods) are associated with the 
generation of each impulse and as a consequence of these 
changes, the firing frequency is modulated by the value 
of suprathreshold excitation. 

vi) in impulse is generated by a cell with some delay (so- 
called synaptic delay) after the arrival of suprathreshold 
stimulus. The extent of this delay depends on the 
intensity of the stimulus and on the state of membrane 
excitability at the given moment. 
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Some other properties attributed, to nerve cells include - 
variations in the weights (effectiveness in inducing the 
impulse) of synaptical inputs and arbitrary adaptational 
changes in the threshold. 

Most neuron models can be classified in two primary 
catagories - (i) those modelling patches of the excitable 
membrane and (ii) those representing input-output relationships 
for the whole neuron. 

3.1 Membrane Patch Models 

The earliest physical model of excitable membrane and 

r p 1 

axon is the ionic model using M iron-nitric acid system. L J 
Most modern models are however, electronic in nature. One 
branch of electronic models is centred around the development 

m 

of 'neuristors' based on a concept introduced by Crane. J The 

r i 

major emphasis in this kind of studies L J has been on a 
device having the form of a one dimensional channel along which 
signals may flow, the signals taking the form of propagating 
discharges having the following properties - 

(i) threshold stimulatibility 

(ii) uniform velocity of propagation 

(iii) lossless propagation 

(iv) refractory period following a discharge 
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Most of the other models are based on the mathematical 

r 7 t 

model for a neural membrane given by Hodgkin--Huxley. J Many 
complex mod els ^ ^ have been devised, which do incorporate 
subthreshold passive responses as well as action potential 
generation, often to the point of identifying the transmembrane 
ionic currents involved. 

Most of these models are relatively complex a.nd are 
thus uneconomic for extensive neural networks. 

A relatively * Slimple model incorporating both active 
and passive properties may be constructed' based on the 
conventional representation of an axon, as shown in fig. 3.1. 

It breaks the axon into discrete segments and represents the 
axoplasmic resistance by Ra, the extracellular resistance by 
Re, the transmembrane leakage resistance by R m transmembrane 
capacitance by C m and the active process involved in the 
generation of the action potential by a nonlinear time varying 

T 1 2 1 

resistance . Brockman L J has proposed a simple two-- 
transister circuit simulating R m ,C m and R^. 

3-1.1 Brockman's membrane model 

The basic circuit as proposed by Brockman is shown in 
fig. 3.2, which is essentially a collector coupled controlled 
astable with a large base resistance R£ at . To explain 
the operation of the circuit, let us first consider the 
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condition for d-c steady state to exist and then establishes 

bias such that Q. is in cutoff while R- biases Q 0 in saturation. 

1 o 2 

The impedance looking into terminals A-B is then a capacitance 
(in series with the small saturation resistance of C^) in 
parallel with a resistance given by 

R m = E 1 // < R 2 + < R 3a// E 3b^ < 3 -’> 

As long as the transmembrane voltage V m does not rise 
to a voltage where begins to conduct, the membrane is 
adequately represented by these passive elements. 

let us consider the effect of a gradual increase in V . 

As V increases the threshold 7^ ; given by the equation 

V T = [R 2 Hy ^)] V (H 3aV/ 

- [(R 2 //E 3 b )V 00 tH 2 + (E 5 a /R 3 b ))]/[(R 3 a+ (R 2 //R 3 b ))(R 3 a ^H 3b )] 

( 5 . 2 ) 

begins to conduct and because of the loop gain being 
more than one, a regenerative action takes place. The rise in 
voltage at point A corresponds to the first phase of an action 
potential (depolarisation of the membrane). The duration of 
the pulse at A is determined by . The fall in the voltage 

corresponds to the r epolarisation phase. By a proper choice of 
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and C m ratio, the hyperpolarisation phase also can be 
observed in the form of V undershooting below the resting 
value. 

If the model is excited by a current I . V starts 

ex m 

rising towards I . R with a time constant R C . In case 

ex m mm 

I ex there is a firing the moment V exceeds V^. 

Following the pulse, the process of charging starts again 
and a volley of impulses result. 

Brockman studied the behaviour of an axon simulated by 
connecting sections of this membrane patch model with R^= 10K 
and Rj, = 0 (allowing the use of a single power source) and has 
reported that the phenomena of electrotonic and action 
potential propagation are reasonably simulated by this model. 

This model can be used for simulating myelinated as well 
as unmyelinated fibres. In the case of myelinated fibres, this 
active network is placed at the points corresponding to ±he 
nodes of Ranvier, connected by passive R-C ladders. 

This circuit exhibits relative refractory period and 
thus is capable of frequency modulation. However, it does not 
exhibit an absolute refractory period. The impulse waveform ■ 
is not a good approximation of the action potential and the 
threshold is not accurately adjustable. Further, -the amount 
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of undershoot at the end of the impulse depends on the pulse 
width (determined by ) and hence the refractory period 

and pulse width cannot be varied independently . 

3.1.2 A modified membrane patch model (version - I) 

The simple model proposed by Brockman may be modified 
to incorporate an absolute refractory period, which can be 
achieved by inhibiting the input for certain period immediately 
after the end of the pulse. This has been done by including 
a transister (as shown in fig. 3. 3) in the model. Because 
of its base circuitry, is off under d.c steady state. At 
the rising end of pulse, diode D conducts for a short time, 
quickly discharging capacitance as falls. The rise is 
rounded because of the effect of Q . At the end of pulse, the 
regenerative process is somewhat slowed down because of effects 
of C 2 and- C^, but if R^ 0^ , Rg 0^ < < R m C ffi , it will still 
result in a termination of the pulse with an undershoot which 
is -arrested to ~ 2V. ^ and quickly comes back. The rise in 
turns on which prohibits excitation as the stimulating 
current is shunted by its collector. It comes out of this 
state after a time, corresponding to absolute refractory period, 
which is determined by (R^ + R^)0^. 

In addition to the incorporation of an absolute 
refractory period, the circuit has the additional advantage 
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at th& cad o t absolute refractory period, v m i 0 a ^ ^ CBS * 

irrespective of ratio of 0 an d C. an d thus the puls® width 

ml 

and relative refractory period may, be varied ind ep el1 ^ eI1 ily • 

3.1.3 Modified membrane patch model (version II) 

fhe modified version discussed earlier does incorporate 
an absolute refractory period, but the waveform of ^ iie im - pu ^~ se 
is not a good approximation of action potential. 4n 
approximation of action potential waveform may be achieved by 
using separate sections in the circuit, controlling the vari °us 
phases of the impulse activity. A realisation of such a 
circuit, while still retaining the basic properties earlier 
model has been developed as shewn in fig . 3-4. In the d * c 
steady state Q 1 , and are off and Q 2 and Q 4 are on. Thus 
for subthreshold signals, the properties are the same as 
discussed for earlier two models. ~Q 2 and Q^~Q 4 f° rm 
monostables. As v m exceeds the threshold due to a 
suprathreshold current excitation, the first monostabl e 
diggers and voltage at A goes high. Is the rise in ^ c2 
(slowed by R ^ C 2 ) reaches a threshold the second monostatle 
trigger s, which after a delay (due to effect of ) furns 
°n and thus process of depolarisation starts. remains on 
for the second monostable pulse duration thus providi n & an 
absolute refractory period (determined by Rg 4 0 1 ). 
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Except for the later portion of falling edge, the 
impulse -waveform is a good approximation of action potential. 
The falling edge may he made more smooth by introducing an 
R-C ladder at the base of Q,- . This circuit also suffers from 
the lack of an accurate adjustment of threshold. 


3.2 Neural Cell Models 

When one is primarily interested in the function of a 
nerve cell in information processing, rather than in the 
mechanisms underlying this function, many of the complexities 
related to active processes and shape of waveform of action 
potential may be ignored while modelling. 

ri3i 

Most of the models J simulating the input-output 
relationship for the whole neuron rather than small areas of 
membrane separate the input terminal from the point where 
modelled transmembrane voltage is developed. The earlier 
efforts in this kind of modelling were making use of discrete 
components. The basic problems associated with these models 
are the complexity of the circuit and a large number of 
undetermined parameters. In the more recent models- 
integrated circuits have been used. However, the no. of chips 
(typically 3 to 10) and the interactive nature of parameters 
are the basic hurdles in using them for extensive neural 


studies. 
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The emphasis in the present work has been in constructing 
■a simple model with a reasonable functional representation. 

For the purpose of simulation, the various processes associated 
with neural activity have been divided in three sections, and 
a circuit simulating each has b een d eveloped (fig. 3. 5). The 
first circuit has many inputs and simulates the properties of 
temporal and spatial summation as well as adaptation. Its 
output is graded and hence is termed as generator potential 
corresponding to a receptor cell output and the effect of 
synapses. The second circuit is the action potential 
generation block which converts the graded generator potential 
into a volley of impulses. The third circuit is an active 
delay line simulating the lossless propagation along the axon 
length. 

3.2.1 Generator Potential Block 

For the simulation of various analog processes viz ~ 
adaptation, temporal and spatial summation, circuit shown in 
fig. 3.6 has been developed. Circuit in fig. 3.6a simulates 
the effects due to an excitatory input channel. Combination 
of R^ , R^ an d simulates the process of adaptation. The 
process of temporal summation (integration) is simulated by 
the combination of R^, R^, C^. The build-up time constant 
corresponding to the release of the chemical transmitter at 
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Pig. 3*6 Generator Potential Block 

(a) Circuit oh owing one cxcitating input 

(b) JiTect of eheptation 

(c) Circuit. showing' one cxcitating read one inhibitory input. 



Pig. 3*7 Scheme for action potential generator block. 
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the synapse is controlled by R^ while the deaay time 
constant, corresponding to the decomposition of chemical 
transmitter in the synaptic gap is controlled by R^ C^. Fig. 
3.6c includes the effect due to an inhibitory input channel 
also. The components and the time constants are in 
correspondence with the ones for excitatory inputs. A number 
of excitatory and inhibitory inputs may be connected together 
with 0^ as the capacitor on which temporal and spatial 
summation of charges take place. The transistors interfacing 
the adaptation and summation may he looked upon as synapses with 
unidirectional transmission property. 

3.2.2 Action potential block 

Action potential generation can be simulated by a simple 
scheme shown in fig. 3.7. The input voltage (generator 
potential, the output of the first circuit) is continuously 
compared with a threshold V^. When exceeds comparator 

output 7 ] goes high and pulse is generated at the monostable 
output. Falling edge of pulse triggers the threshold circuit 
and its output goes high ; and hence comparator output goes low. 
The threshold remains at high level for a certain duration 
corresponding to absolute refractory period and then starts 
decaying towards its steady state value with certain time, 
constant. Another pulse is generated as soon as has fallen 
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below V^. This process results in an encoding of the 
magnitude in the form of frequency of pulses. 

The threshold circuit may be realised by a negative 

edge triggered monostable determining the absolute refractory 

period with a circuit giving a decaying threshold after this. 

period as shown in fig. 3.8. Steady state threshold is 

determined by R and R while decay rate is controlled by 
* x y 

C (R x /7 R y ). 

This scheme may be realised in a number of ways using 
commercially available monostable and comparator ICs. However, 
it is possible to use timers in place of the monostables as the 
inputs to them remain in steady state for the duration of 
quasi-stable state. As the time-durations involved are of the 
order of a ms, CMOS inverting gates have been used. The 
comparator also may be made with four CMOS inverting gates 
[Appendix], Thus the whole scheme has been realised with one 
CMOS chip with six inverting gates and one PHP transistor as 
shown in fig. 3.9. 

Steady state threshold V^ r is given as 

V Tr ~ / S 5 + V T cc 

(neglecting the effect due to Rg) 
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Pulse widths of the two timers, T^ and are gi ven 


as 


T 1 = Vl X n < V 0H /V ^ 0-7R 1 C 1 
T P = V 2 l n (V QH /V t ) - 0.7 R 0 C. 


2 2 


(3.4) 

(3.5) 


where V = threshold of the CMOS gate. 

Por the absolute refractory period, the threshold repairs 
^cc ” ^"cES as ^ is driven "5° saturation. To ensure t uat Q 
goes to saturation during this period, 


at 


^ f> h 


(5-6) 


and 


where p = dc current gain of Q. 

Following the absolute refractory period, V 2 goes high ; 
hence Q cuts off. decays to its steady state value with 

time constant t. 


T 884 (®3 V 0j 

Threshold voltage is given as 


(3.7) 


V T = ((7 


cc 


)- e 


CES ' Tr 


t/r 


+ 7, 


Tr 


(3.8) 


-Another pulse is generated after a time T^, when 7^ just 


falls below 7^. Thus we have 


T, 


R 3 C 3 l n ( (V nn - 7 rp J/(7 n _ - 7^)) 


00 v Tr' /vv G Tr ' 


(3.9) 
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The impulse frequency is given as 

f = 1 / ( t 1 + T 2 + t 3 ) (3. id) 

Eor Vq -• V^, no impulses are generated, while as Vq approaches 
V CC ~ V CES kig&est value of threshold), the impulse 

frequency reaches a maximum f ~ 

f ma* = ’/< T 1 + V 

Pig. 3-10 gives an impulse frequency Vg plot. It is 
to he noted that the maximum threshold is finite (V^ - Vp-gg) 
and for any input above this value impulse generation will 
not take place as the comparator output will always remain high. 

3.2.3 Active delay line 

The axon may be simulated by an active delay line on 
which impulses propagate without attenuation. This may be 
achieved by- several timers (using CMOS inverting gate) in 
cascade or sections of RC ladders with timers at the points 
corresponding to the nodes of Ranvier. Each timer may be used 
for a maximum delay equal to the pulse width. 



CHAPTER 17 


THE NERVOUS SYSTEM 

In order to survive, the living organisms must adopt 
themselves to changing conditions. Higher organisms have 
developed a centrally controlled system which creates 
conditions favourable to the organism as a whole. There are 
usually several centres, with hierarchical relations among 
them. According to the manner in which signals are 
transmitted, two typical systems may be distinguished. 

(a) The glandular system, where signals are 
transmitted through blood circulation system by chemical 
compounds. 

(b) The nervous system, where signals are transmitted 
over nerve fibres by electrical pulses. 

In addition to these systems, there are automatic 
regulation systems of individual organs. All these systems 
are interactive and a close cooperation between them forms 
the basis for survival. In this work, the interest is 
focused on the nervous system. 

The nervous system is the supreme control organ and 
is found only in living organisms at higher stages of 
evolution. It processes an enormous amount of information 
obtained from the environment as well as from internal 
organs. In vertebrates, . it is made up of the brain, the 
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spinal chord and the many nerve processes that pass between 
these two structures, as well as muscles, glands or receptors 
which they innervate. 

There are several schemes of classifying the nervous 
system, the three possible bar~its for classification being 
the anatomy of the system, its function and an appropriate 
combination of the two. In one of these schemes the brain 
and spinal chord together are termed as the central nervous 
system (CNS). The nerve processes outside it form the 
peripheral nervous system. 

4.1 The Central Nervous System (CITS) 

This is the supreme coordinating nervous system. It 
maintains the functional integrity of higher living organisms 
and controls their behaviour so that they are adapted to the 
given conditions in the most advantageous manner. 

The simplest activity mediated by CITS is an 

elementary reflex which involves a receptor (sensory organ), 

afferent nerve fibres (transmitting sensory information to 

CNS), certain parts of CHS, efferent nerve fibres 

(transmitting control information from CITS) and an effector 

(executive or control organ). Reflexes can be classified 

(" p 3 1 

from several points of views. * J According to whether 
one is concerned with a contact with the environment of the 
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organism or its internal organs, reflexes are divided into 
somatic and autonomic ones. According to the part of CITS 
involved, they are divided into spinal and eerebral reflexes. 
According to whether they are acquired or innate, they are 
called conditioned or unconditioned reflexes. All the 
activities of CITS may he looked upon as a set of different 
kinds of reflexes, interacting with each other in a complex 
manner. 

r3-5l 

CITS can he roughly divided into seven basic parts J 

(i) spinal chord 

(ii) medula oblongatta (spinal bulb) 

( 111 ) cerebellum (little brain) 

(iv ) mesen cephalon (mid brain) 

(v) diecaphalon (between brain) 

(vi) carpus striatum 

(vil ). cerebral cortex 

These parts are in the order of both morphological 
and functional hierarchy. 

Some lower control functions are directly performed 
by the spinal chord, whereas more complicated functions are 
left to the higher parts of dTS. In the spinal chord, all 
the reflexes without any exception, are ■unconditioned. 



The spinal bulb performs some important functions - 
control of the action of blood circulation system, the 
respiratory control, the digestive system, defense reflexes, 
movement of facial muscles. It also plays important role on 
normal body posture. The little brain controls the 
coordination of more complicated movements. The midbrain 
mainly acts as an interface and in addition to it controls 
eye movements and mediates unconditioned auditory reflexes. 
The between brain has two parts - (a) Thalamus and (b) 
Hypothalamus. Thalamus processes all excitations resulting 
in various sensations (cerebral cortex also participates in 
them). Hypothalamus coordinates vegetative functions -• blood 
circulation, temperature regulation, respiration etc. Corpus 
Striatum serves as a means of transformation, which controls 
the goal seeking function of skeletal muscles according to 
the orders from the cerebral cortex: Cerebral cortex performs 
all higher nervous activity. The fundamental element of 
cerebral cortex is conditioned reflex. 

I 

4.2 The Peripheral Nervous System 

The nerve processes in the peripheral nervous system 
may be classified into afferent and efferent divisions. Cell 
bodies of afferent neurons are outside CHS. Regardless of 
whether the initial stimulus is picked up by exteroceptors 
(external sensors) or interoceptors (signalling internal 
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conditions ), the action potentials travel along neurons 
which are structurally similar. These are also referred to 
as sensory neurons. 

The efferent division may he sub-divided into somatic 
and autonomic nervous systems. The somatic nervous system 
innervates skeletal muscles and the autonomic system 
innervates smooth and cardiac muscles and glands. 

The cell bodies of the neurons in the somatic division 
are located in groups in CNS and the axons pass directly 
without any synapses to skeletal muscle. Their excitation 
always leads to a contraction of muscle cells, there being 
no inhibitory ones. These cells are also called motoneurons; 

The fibres of the autonomic nervous system have 
synaptic interconnections in the intermediate region between 
CNS and the neuroeff ector junctions. The autonomic nervous 
system can be further subdivided into sympathetic and parav 
sympathetic components. They leave CNS at different 
levels and for most part are functionally reciprocal. The 
sympathetic nervous system tends to mobilise the body for 
emergencies, whereas the parasympathetic nervous system tends 
to conserve and store bodily resources. 

4.3 Information Processing, Learning and Memory 

Neurons encode the varying magnitudes of stimulus 
by the process of frequency modulation. This has a natural 



40 


immunity to noise. Although, the neuron is digital in action 

(action 'potential being all or none phenomenon), nature has 

discarded the much more efficient and compact positional 

numbejr system in favour of the mere primitive method of 

counting. The basic advantage of this is in the greater 

safety against errors and greater simplicity of encoding - 

[7] 

decoding. It also results in a large amount of redundancy. 

Sensory organs supply CHS with a huge amount of 

Q 

information (according to some estimates - 10 bps) about the 
environmait and internal state of the organism. It includes 
not only important, but also redundant information. CNS can 
not record and process such an amount of data. The sensory 
channels themselves encode these data, reducing the redundancy 
by filtering the parameters suitable for furtheh processing. 

The processing inside OH’S takes place at all the hierarchical 
levels. 

Learning is a kind of plasticity of the nervous system. 
The nervous system, specially in mammals, can store information,, 
can modify response to stimuli and can recover from certain 
irreversible structural damages. The processes of memory, 
logic and arithmatic involved are not well understood. 

Several theories relating learning to structural changes at 
the synapse or the chemical alterations within the nerve 
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cells have been proposed. However, the problem of total 
organisation of nervous system for memory and learning is 
unsolved. Theories of memory involving a physical change 
as the memory trace are of two sorts. Short term memory is 
explained in terms of active electrical processes, such as a 
reverberating circuit - a closed chain of neurons reactivating 
one another. Long term memory is thought to be the result 
of a more permanent alteration in brain chemistry or synaptic 
structure. ^ 



CHAPTER V 


THE CARDIOVASCULAR SYSTEM 

The activity of the organs of cardiovascular system- 
the heart and the blood vessels , ensures a continuous flow 
of blood in the organism. Blood circulating in the organism 
performs a number of functions - gas exchange (supplying 
oxygen and carrying off carbon dioxide), nutrition of tissues, 
control of body temperature etc. 

The flow of blood in the vessels is due to pumping 
action of the heart. As the requirements of blood for the 
entire body or its parts vary continuously, the heart and 
blood vessels are subject to continuous regulation in closed 
loops interconnected via the nerve centres. 

5.1 The Heart 

A simplified diagram of the human heart from an 
engineering point of view is given in fig. 5.1^^ The right 
half of the heart pumps blood from the body into the lungs 
and the left half from the lungs to the body. These two 
subdivisions are called pulmonary or lesser and systemic or 
greater circulation respectively. 

The structure of the two parts is the same - each 
consists of an atrium and a ventricle, separated by a valVe 
which passes blood only from the atrium to the ventricle. 
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Ft -:. -5.1 A simplified digram of the human heart 



Fig. 5.2 Schematic representation of changes of piaseure 
on (A) right and (B) left heat associated with 
the different phases of cardiac cycle (at 75 hprt.) 
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The outlet from the ventricle is also closed by a valve 
permitting only an outward flow of blood. 

The rhythmic pumping of blood is performed by 

alternate rhythmic contraction (systole) and relaxation 

(diastole) of. the muscular fibres that form the walls of the 

atria and ventricles. In normal physiological conditions 

systole and diastole occur in a definite coordination and 

r p i 

constitute the cardiac cycle. 1 J fig. 5.2 shows the various 
phases of the cycle alongwith pressures , which will be 
discussed later. 

/ 

Consecutive contractions of the heart follow each 
other because of the automatic mechanism of its own. The role 
of the nervous system is only in weakening or strengthening 
the heart action. 

The automatic action is made possible by centres 
where self excitation occurs. These are the areas of 
muscular tissue of a special structure called nodes. The 
stimulation arises initially in the Sino-atria.l (S A) node 
located in the front part of the right atrium near the 
arifices of vaaaa cavac. This spreads through the muscles 
of atria in about 100 mS and reaches the atrio-ventricular 
(AV) node, the sole muscular connection between the atria 
and the' ventricles. Upon passage through this node, 
excitation travels through specialised right and left bundles 
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and then through the arborised PurkinJ e fibres in about 
100 mS. This system serves to initiate electrical activity 
in the v entricular musculature* The two ventricles contract 
simultaneously and the contraction persists for about 300 mS, 
normally followed by a relaxation period of 500 mS (with 75 
beats per minute). The electrical conduction system is shown 
in fig. 5.3 J 3 -* 

The heart rate is basically determined by the 

activity in the SA node which operates, in a sense, as a free 

% 

running multivibrator. The rate, however, is modified by the 
competing effects of sympathetic and parasympathetic nerves 
coming from CFS. The SA node is the main pacemaker of the 
heart. In the case of a failure of this main pacemaker, the 
AV node takes over the role of pacemaker and hence is known 
as the latent pacemaker. The contractions of the atria and 
ventricles may not then occur in their usual sequence, but 
almost simultaneously because stimulus from the AV node 
reaches the atrial and ventricular muscles almost at the same 
time. Ventricular contraction persists even aftej» the A-V 
node has been separated from .tho lower sections of the 
conduction system by ligation. The function of pacemaker 
is then performed by the Purkinje fibres lying in the right 
and left ventricles. 
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The frequency of automatically originating impulses 
shows the degree of a pacemaker's activity. The rate of 
cardiac contraction in a normal man at rest averages 70 — 75 
beats per minute (bpm). It changes with the physical 
activity due to actions of sympathetic and parasympathetic 
nerves and ranges from 60 to 85. With AV pacemaker, the rate 
is 40-50 bpm, while a heart working under the lower lying 
pacemakers contracts at an even slower rate (20-50 bpm).^^ 

5.2 Disorders in Conduction and Heart Block 

For proper functioning of the heart, adequate 
correlation between atrial and ventricular contractions is 
essential. This gets disturbed if the conduction of 

r s 6 1 

excitation is impaired. 1 J This is referred to as a heart- 
block and may be due to congenital defects or accidental 
tissue damage. 

An SA block is said to occur when SA node fails to 
initiate an impulse. When the conduction between the SA node 
ventricles is affected, it is called an AV block, which can 
be partial or complete. In a partial AV block, conduction 
from SA node to ventricles takes longer time or some impulses 
fail to reach the ventricles, while in a complete block, SA 
node impulses do not reach ventricles at all. Depending on 
the region in the conduction path (P,Q,H,S in fig. 5.3), 
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different kinds of disorders may be observed. Whan the 
section P is impaired, atria contract under the control of SA 
node, while ventricles contract under the control of AV node. 
There is a partial or complete dissociation between them 
depending on the degree of impairment. In case of failure 
of AV node itself or in the section Q, ventricles run under 
their own low-level pacemakers. Though the beat rate is low, 
the two ventricles remain in synchronism with each other. In 
another kind of disorder, one or both the branch bundles, 
(Section R and S) may be impaired. Impairment of both bundles 
results in the dissociation of the contractions of the atria, 
left ventricle and right ventricle. 

The normal rate of heart rate varies from 60 to 85 
bpm. However,' 1 sometimes a marked variation from this is 
observed. A lower beat rate (40-60 bpm) is known as brady- 
cardia and a higher beat rate exceeding 90 or 100 and some 
times reaching even 150 bpm is called tachycardia. Bradycardia 
is often observed in atheletes at rest, while tachycardia 
occurs during strenuous muscular work and in emotional shocks. 
A disorder leading to an irregularity of the rhythms is 
known as arrythmia. 

Another serious kind of disorder may cause rapid and 
asynchronous contractions of the artrial or ventricular 
muscle fibres, reaching a rate of 400 bpm (flutter) or 
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600 "bpm (f ibrillation ) . With such fast asynchronous 
contractions tho atria and ventricles cannot fulfil their 
function. Atrial fibrillation can persist for a long time without 
danger to life. Tho rhythm of ventricular contractions in this 
case is irregular because only a small fraction of the artrial 
impulses is passed from the atria to them.. Ventricular 
fibrillation is fatal if immediate measures to stop this condition 
are not applied. 

5.3 Blood Circulation 

The actual physiological system for the heart and blood 
circulation system is da own in fig. 5.4. An equivalent 
diagram from the engineering point of view is depicted in 
fig. 5.5 

The atria play the role of a reservoir that collects 
blood flouring from the veins during a ventricular systole. 

The blood passes from the atria to the ventricles during the 
ventricular diastole. The ventricles act as a pump that forces 
blood into the arterial system at relatively high pressure. 

The blood pressure in the atria, ventricle and artries under 
normal condition are shown in fig. 5.2. £8j 

During atrial systole, the muscular constrictions at 
the entrance from veins contract so that blood flows from the 
atria only in the direction of the ventricles. As the 
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Fig ♦ 5.? Circulatory system - engineering point of view 
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ventricles are relaxed during the atrial systole and pressure 
within them is lower than that in the contracting atria, 
blood enters them from the atria. One way passage of blood 
from the ventricles into the main arteries is due to valves. 
The opening and closing of valves is conditioned by the blood 
pressures on the two sides. 

The movement of blood in heart cavities and along the 
entire cardiovascular system is conditioned by pressure 
gradient and the r esistance offered by pathways. The 
greater pressure developed in the left ventricles compared 
with the right is due to more powerful musc3.es of the former, 
which is associated with the fact that the left ventricle 
has to overcome a greater resistance to blood flow in the 
vessels of systemic circulation. The fluctuations of 
pressures in the aorta and pulmonary artery during the 
discharge of blood from the ventricles follow the changes of 
pressure in the corresponding ventricle. The form of the 
arterial pressure pulse changes as it passes through the 
arteries. The walls of the arteries cause damping and 
reflections. The flow of blood from the ventricles is in 
the form of pulse discharges. However, as the arteries 
branch out into smaller arteries with smaller cross-sectional 
area, the pressures and volumes change and the nature of flow 
changes. As the blood flows into the arterioles, the 
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pressure decreases and loses its pulsating character, iis the 
blood returns to veinous system after travelling through 
capillaries, the pressure is down to about 15 mm Hg and 
finally in the veins returning to heart, the pressure is only 
about 2 mm Hg. 

The systolic (maximum) and diastolic (lowest) pressures 
for clinical purposes are measured in the brachial artery in 
the arm.. In a normal adult, the values are typically 120 and 
80 mm Hg respectively. However, they are subject to 
variation with age, climate and other factors. 

The quantity of blood discharged by the ventricle per 
minute is known as the cardiac output. The average value in 
a man at rest is 4.5 to 5 litres, corresponding to a complete 
turnover every minute. Cardiac output depends both on the 
stroke volume (blood discharged from ventricle with every 
beat) and the beat rate. 

The minimum time required for the blood to complete 
both the systemic and pulmonary circulations is known as the 
total circulation time. The circulation time in humans is 
about 27 heart systoles. The actual time taken by the various 
possible circuits differ widely. It is due to the fact that 
the rate of flow along the walls of the blood vessel is lower 
than that along its axis and that extension of all the 
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vascular regions is not the same. The pulmonary circulation 

r 9 1 

accounts for one fifth of the total circulation time. 

A normal blood supply to the organs and tissues 
requires a definite relation between the volume of the 
circulating blood and the total capacity of entire vascular 
system. This is achieved by several humoral and neural 
mechanisms. About 45 to 50 percent of all blood is 
accumulated in the so called blood reservoirs (the spleen, 
the liver, subcutaneous vascular networks and the lungs). 

These reservoirs play an important role in the regulation of 
blood flow. A redistribution of circulating blood occurs 
when an organ is working. Blood supply to the working organs 
is increased at the expense of decrease of the flow to other 
ones. The total volume of blood is maintained due to the 
continuous formation of blood cells and their destruction. 

5.4 Control of Cardiac Activity 

The action of the heart, and the rate and force of 
its contraction vary with the activity of the organism and 
the various conditions in which it finds itself. ITeural 
control is accomplished by impulses sent from CITS along the 
sympathetic and parasympathetic nerves.^ ^ They innervate 

the SA node, the muscle fibres of atria and the AV node. 
Impulses on the sympathetic nerve cause a more rapid 



54 


spontaneous depolarisation of the pacemaker cells during a 
diastole, which leads to the acceleration of heart 
contractions. The action due to parasympathetic nerves is 
the reverse. The nerve centres giving rise to these two 
nerves are in continuous excitation, called central tone. 

In the resting state, parasympathetic elenents are dominant. 

Reflex control of cardiac activity involves many 
nerve centres in CHS. The blood pressure is monitered at a 
large number of points through presso receptors which 
stimulate the different fibres. A higher arterial blood 
pressure causes increased impulse rate on them. This 
increases the tone of the parasympathetic centres, which 
makes the heart beat more slowly. A rise in the blood 
pressure in venae cava o gives rise to a reflex decrease in the 
tone of parasympathetic centres and to the stimulation of the 
sympathetic nervous system as a result of which heart beats 
faster and pumps more blood in the arteries from the veins 
and the pressure in the venae cavae falls. A difficiency of 
oxygen and other chemicals during physical exertion leads to 
stimulation of sympathetic centres which leads to increase 
in cardiac activity. 

In addition to nervous control, cardiac activities 
are influenced by humoral controls and many other factors 


such as temperature etc 



CHAPTER VI 


SIMULATION OP THE CARDIOVASCULAR SYSTEM 

Simulation of the cardiovascular system consists of 
the simulation of two distinct mechanians (i) conduction of 
excitation in heart and (ii) hlood circulation and various 
related processes. 

6.1 A Simple Simulation of the Conduction System 

The simplest simulation of the cardiac output and its 
effect is just a charge pump as drown in fig. 6.1. As the 
activities of different chambers of the heart are not 
distinguished in this model, this allows the simulation of 
conduction of electrical excitation in terms of three distinct 
pacemaker activities in the heart, as follows 

(i) A controlled oscillator representing the SA node - The 
autonomic nervous system (AHS) controls the running frequency 
by two control inputs sent through the sympathetic and 
parasympathetic nerves, 

(ii) An oscillator representing the AV node - The free running 
frequency is lower than that of the SA node. Under normal 
conditions, pulses from the SA node synchronise the AV node 
oscillator to the same frequency as that of the SA node. 
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Pig. 7.1. A simple simulation of cardiac activities 
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(iii) An oscillator representing low level pacemakers in 
ventricles. In the event of an AT block, it runs at a low 
frequency as a free running oscillator. Under normal conditions 
of AT conduction, this oscillator is synchronised with the 
higher pulse rate of the AT node oscillator. 

The output of the ventricular oscillator drives the 
charge pump and there is a release of certain charge on the 
capacitor C corresponding to the heart output with each beat. 

A leakage resistance R^. across the capacitor simulates the 
diffusion of blood. The voltage across the capacitor may be 
taken as a measure of cardiac activity and fed to CUS through 
a sensory neuron. 

6.2 Simulation of the Circulatory System 

The simple model of fig. 6.1, exhibits the three levels 
of automation in the excitation system. But in order that any 
simulation be really meaningful, it must include some kind of 
simulation of the circulatory system as well. An electrical 
simulation can be developed using the following electrical 
analogues for hydraulic parameters - 

(a) pressure - voltage 

(b) flow - current 

(c) volume - charge 
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With these analogues, the friction in blood vessels 
may be modelled by resistance; the capacity of a part by 
capacitance and compliance by inductance. In the heart, opening 
of a oneway valve is decided by the differential pressure across 
the valve, which may thus be represented by a diode. 

For the sake of a meaningful yet reasonably simple 
simulation, the following assumptions and approximations have 
been made - 

(i) The circulating fluid (blood) has no compliance. Thus 
no inductors are included in the simulation, 

(ii) The entire region of any of the heart chambers contracts 
or expands simultaneously (in contrast to the actual case of 
propagation of contraction wave). 

(iii) The stroke volume is constant (in contrast to the 
variations in stroke volume in accordance with the needs of 
the body. 

(iv) The distributed circulatory system can be represented 
by a finite no. of sections. 

A simulation of the circulatory system has been developed 
as shown in fig. 6.2. The contractions of various chambers 
are simulated by pulse drivers with outputs in correspondence 
with the actual forces of contraction. The muscular constriction 
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Pip. 6.2 Simulation of circulation system 
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at the inputs to the atria are represented by diodes shunted 
by resistors. The pulmonary circulation is represented by two 
sections of R-C ladders, while systemic circulation is 
represented by five such sections. 

Before the last section of systemic circulation a 
regulating reservoir has been introduced. Its function is to 
regulate the amount of charge (blood) in circulation. It 
incorporates the functions of blood reservoirs and the organs 
taking part in continuous formation and destruction of blood 
cells* 

Effects of bleeding as well asctransfusion can be 
simulated by leaking the charge or introducing the charge 
somewhere in the circulatory system. 

6.3 A Detailed Simulation of Conduction System 

In view of the model of the circulatory system proposed 
in the preceeding section, certain features have to be 
incorporated in the simulation of conduction system so as to 
make them compatible with each other. Eor the sake of avoiding 
complexity, the various pacemakers are treated as localised 
points (which is justified because in the simulation of 
circulation, propagation of the contraction wave in the chambers 
of heart has not been considered). Further the effects of , 
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sympathetic and parasympathetic nerves in the autonomic nervous 
system is assumed to have effects on SA node only. Thus the 
simulation has to incorporate the following features - 

(i) Separate pacemaker activities of the two ventricles, 
(ii) Bidirectional paths envolving certain delays Between 
the various nodes, which would allow proper conduction among 
the contractions of the various chambers. 

For simulating the various node oscillators as well as 
delay paths between them the circuit shown in fig* 3.3, proposed 
for simulating the membrane patch has been used. This relatively 
simple circuit features an absolute refractory period and can 
be made self oscillatory by increasing the resting value of 
above the threshold, which can be readily accomplished by 
connecting a resistor from the supply to the input. By changing 
this bias, different free running frequencies, • corresponding to 
different nodes may be achieved. Eor the bidirectional delay 
paths between them, same circuit with no self excitation has 
been used. The whole conduction system isthus simulated as 
shown in fig. 6.3. 

The excitatory and inhibitory controls on the SA node 
activity have been simulated by adding two transistors Q 1 and 
Q 2 to the basic circuit. Each impulse on the sympathetic nerve 
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For simulating the different forces with which the 
muscles of different chambers contract, the pulse heights to 
the inputs of the capacitors (representing the heart chambers) 
have been made adjustable. As the two atria do not have separate 
pacemakers, normally there is no conduction disorder and there 
is not much of difference of contraction forces of the muscles, 
a single output from the point before the delay between SA and 
AV nodes, has been taken to represent the muscle forces of both 
the atria. The two ventricles have distinct paths of conduction 
from the AV node (the left and right His bundles), different 
forces of contraction and separate low level pacemakers. All 
these features are incorporated in the simulation. 

6.5 Feedback to CHS 

The two kinds of parameters related to the cardiovascular 
activities, which are sensed and • fedback to CHS are - 
(i) Pressure levels at different points in the 
circulation system, 

(ii) levels of oxygen as well as various chemicals 
required for metabolism in the tissues. 

As it is not possible to include anything analogous to 
the level of oxygen and various chamicals in the charge, the 
analogue for blood, we assume that these levels are constant, 
thus ignoring the effects of variations in breathing and 
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various humoral processes. The reserve of oxygen and other 
vital chemicals is modelled by the charge on the capacitor C in 
a charge pump as shown in fig. 6.4. The £L ow rate is sensed by 
sensing the voltage across a resistor in the circulation path. 
In order to avoid Op Amps, the difference amplifier has been 
made ua hg CMOS inverting gates (vide-Appendix) with transistor 
buffers at the inputs, although it does have a finite leakage 
current leading to drift, A charge proportional to the flow 
rate is released on the capacitor C, the voltage across C 
representing the level of the reserves and this can be fed to 
CHS through a sensory neuron. 

The effects of physical exertion, resulting in depletion 

\ 

of the reserves is simulated by a variable leakage resistance 
across 0, 
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CHAPTER VII 


SIMUIuiTICH OP THE HER7CUS SYSTEM 

The processes involved in the peripheral nervous system 
may he simulated hy using receptor models and a suitable 
interconnection of neural cell models. Simple unconditioned 
reflex involving spinal chord can also be simulated -with one 
or more synaptic junctions in the spinal chord. 

It is possible to achieve the processes of logical 

Vj 

functions and short time memory by using simple neuron models. 1 * 
However, the approach involves a tremendously large no. of 
neuron cells to exhibit any kind of conditioned reflex activity 
or learning process. Therefore, in the present work higher 
activities of CHS have not been simulated by using a neural net 
but- a microprocessor based system has been used instead, Phis 
allows one to use the same hardware for simulating various 
kinds of conditioned reflexes and learning processes simply 
by making appropriate changes in the software. 

Por the peripheral nervous system the neural cell model 
discussed in chapter 3 has been used as a basic module. 

7.1 The Brain Simulation 

The brain simulation has been realised with a micro- 
processor system with several input ports (simulating afferent 
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division) and output ports (simulating efferent division)* 

The processor is connected to a RAM, which serves as a memory 
during processing and a ROM, which stores the program for 
simulating the various processes. 

Because of the processing time involved and the hardware 
cons iderat ions, the no. of afferent and efferent nerves was 
restricted to six and four respectively. 

A general picture of the organisation is dhown in 
fig. 7. 1 . 

The processor handles the ports and does the processing 
in real time with the help of a timer. The timer keeps 
interrupting the processor at certain regular intervals 
(cycle-time). After each interrupt, the processor reads the 
data from' the input ports, does the processing, serves the output 
ports and then goes to a wait loop, a similar cycle is initiated 
hy the nert interrupt. 

Due to the availability of various software development 
facilities and support chips, the processor chosen is Intel 
8085. It is operated at 2.5 MHz clock (provided by on the 
chip clock generator and an external 5MHz crystal). 

7.1.1. Input int erf ace 

Since Intel 8085 is a parallel 8-bit processor, the 
impulse rates on the afferent channels have to be converted 
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to an 8-bit parallel data word. This has been achieved by 
giving the impulses on the channel as a clock to an 8-bit 
counter. The counter output is read by the processor as the 
strength of the signal an the channel and the counter is cleared 
for accumulating the counts in the next time slot. 

7.1.2 Output interface 

The output on each efferent channel also has to be in 
the fjrm of of a series of impulses; the frequency being 
proportional to the signal strength. This has been achieved 
by using an 8-bit pulse rate multiplier (PBM). The output data 
is latched onto the PBM input and a clock derived from the 
system clock is given as the PBM clock. 

7.1.3 Cycle time 

In choosing the cycle time, one has to make a compromise 
between the temporal resolution of sensory signals and the 
complexity of the learning processes that can be simulated. 

A value of 100 ms has been chosen, as a value higher than this 
may result in unacceptable loss in time resolution. However, 
this interval is decided by a software controlled timer (on 
Intel 8156 chip which also provides 256 bytes of RAM and 2 output 
ports) and one has the option of choosing his own cycle time 
by modifying the setting of the timer. 
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7* 1.4 Memory & I/O map 

A direct I/O scheme has "been chosen. A circuit diagram 
of the system has been shown in fig. 7.2. The memory and I/O 
map is given in Table 7.1. 

TABLE 7.1 


(A) Memory map - (iO/M = 0) 


Address (Hex) 

Selected Device 

Remarks 


0000 - 07iT 

2716 

(EPROM) 

ADq rj is latched 

using 




ALE to “provide 

A 0~7 for 




2716. D 0 __ ? of 2716 is 
connected directly to 




AD 0-? of 8085 


0800 - 08PP 

8156 

(bam) 

AD q 7 and all controls 




of 8156 are connected to 
the corresponding lines 




of 8085 


(B) 10 map - (IO/M 

= 1) 




Port Address (Hex) 



Port Assignment 


0 



Input port 1 

(si) 

1 



Input port 2 

(S2) 

2 



Input port 3 

(S3) 

3 



Input port 4 

(S4) 
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4 

5 

6 

7 

8 

9 

A 

B 

C 

D 


Input port 5 (S5) 

Input port 6 (S6) 

Output portl (Ml) 

Output port 2 (M2) 

COR of 8156 

Output port3 (01 )-PA of 8156 
Output port4(C2)_PB of 8156 
Unused 

lower timer Byte of 8156 
Higher timer Byte of 8156 


7.1.5 Program structure 

The program of the software to Be loaded in 2 K ROM 
should have following structure - 

(1) INITIALISATION - The timer and I/O ports of 8156 are 

\ 

programmed. Various registers and RAM locations are cleared 
or set as desired By the programs. 

(2) WAIT - A wait loop. The processor .remains in the 

mwBA. 

loop nifti*— I a hardware int errupt on RST 7.5 (sent By the 
timer) is received and goes to the processing routine 
RS75. 

(3) RS 75 - Service the I/O ports and do the processing. After 
completion of the processing processor goes Back to the 
wait loop in (2). The total time taken By the routine 
must Be less than the cycle time. 
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fig. 7.2 c Microprocessor based system for f feraln* sinmUrtlift. 
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(4) TR AP - Whenever unmaskable hardware interrupt TRIP is 
received, it indicates a catastrophic sl. tuation and the 
processor immediately starts executing this routine. 
Processor may return to nomal mode after TRAP line 
goes low. 

7.2. Peripheral Nervous System Simulation 

Por simulating the afferent division of peripheral 
nervous system, the model developed for a nerve cell in section 
3.2 has been used. However, the axonal delay simulation has 
been omitted because for simulating any reasonable amount of 
delay, a large number of timers are needed and it is really of 
not much consequence from the point of view of simulation. 

The first two sensory channels have one excitatory input 
each. ' The third channel has one excitatory and one inhibitory 
inputs, while fourth channel has one inhibitory and two 
excitatory inputs. 

The remaining two channels do not have the generator 
potential block, i.e. analogue processes, viz - adaptation, 
temporal and spatial summation etc. are not incorporated. These 
are suited for sensing the internal parameters, viz - blood 
pressure and the level of reserves of vital chemicals. An 
unconditioned spinal reflex has been incorporated by a synaptic 


/ 
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function after the first cell (corresponding to a synapse in 
the spinal chord) which excites a motoneuron (only the action 
potential block). 

The efferent division of PUS has five members. Two of 
the outputs from the'brain' form the autonomic division and 
control the cardiac activity by innervating the SA node. The 
other two outputs from the 'brain 1 and the motoneuron output 
in the spinal reflex form the somatic division. In actual case 
one motoneuron excites many muscle fibres, the response of 
individual fibres being on all or none principle. The effect 
is the result of the fact that a particular muscle is excited 
by a large no. of nerve fibres. In the present case, the 
interest is in seeing the effect of change in impulse rate on 
a motoneuron, over a muscle output and therefore, in place of 
simulating individual fibres, a circuit which converts the 
impulse rate to a corresponding dc level (shown in fig. 7.3) 
has been used. There are two such muscles. One is excited by 
the spinal reflex as well as one output involving brain, and 
by breaking the paths; effects of spinal reflex may be seen. 

The other one is excited by the output involving brain only. 
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(a) Muscle with one eroitation input 



IL - 100K 

(b) Muscle with two excitation inputs 

C-10 f 


Fig. 7.3 Muscle Simulations 



CHAPTER ¥111 


SYSTEM ORGANISATION AND OPERATION 

The simulations of nerve cells, the cardiovascular system 
and the nervous system as developed and discussed in chapters 
III, VI and VII have "been made to form a system, which may be 
used to perform a variety of experiments. In this chapter an 
overall organisation, printed circuit card details, operating 
details and suggested experiments are presented* 

8.1 Overall Organisation 

The system has been made in three cards - ( i 5 ONS (brain) 
simulation, (ii) cardiovascular system simulation and (iii) the 
peripheral nervous system (PNS) simulation (alongwith muscle 
simulations and one spinal reflex involving a synapse). 

As the simulation of CITS is based on a microprocessor, 

the internal signals are of little interest from the bioelectric 

point of view. Hence for the sake of simplicity, all the data, 

control and address buses have been kept internal to the card. 

The external connections to the card consist of the power supply 

lines, six serial input lines (afferent nerves S1-S6), four 
serial output lines (efferent nerves: M1,M2,C1,C2) and two 

system control lines (RESET and TRAP). 
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The PH'S and cardiovascular system simulation cards are 
both mounted on the back of the front panel* with PC test pins 
coming out through holes in the panel. The various points of 
interest have been made easily accessible to the user through 
these pins. The cardiovascular system simulation card has only 
the power supply connections while the PHS simulation card has 
twelve lines from the CHS card (S1-S6, Ml, M2, Cl, C2,TRAP, 

RESET) also connected to it through a connector. There is a set 
of potentiometers and jacks on the panel which have been provided 
for different purposes. The whole set up is such that all the 
cards and the panel may be easily separated. 

The entire organisation of the system is printed on the 
front panel, shown in fig. 8.1. The internal connections 
between various blocks are shown by solid lines, while the ones 
to be made by the user are suggested by dotted lines. All input 
points use black pins, while output points use red pins. In 
cases where one msy be interested in seeing the signals even 
after making the external patch card connections, additional 
test pins (blue) have been provided -s they are in parallel with 
corresponding output pins in many cases, one should be careful 
in using any of such points as inputs. Grey pins have been used 
for accessing signals representing the blood pressure in the 
circulation system and the forces of contraction of cardiac 
musculature. The monitoring of blood pressures should be done 
using high impedance probes. 
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8.2 PC Card details 

The details of the three PC cards are discussed in the 
following subsections. 

8.2.1 CITS Card 

The schematic of the circuit on this card is given in 
fig. 7.2. The functions of various components are listed in 
table 8.1. 


Table 8.1 
list of functions and nos. 


Components Functions 


of components on CNS card (pc 1 


Nos. 


1) Intel 8085 

2) Intel 8156 

5) Intel 2716 
4) Intel 8282 


8-bit microprocessor 1 

256-byte of RAM, two 'l 

8-bit output ports and 
one programmable timer 
2E byte EPROM 1 

(a) 8-bit latch/tri-state 6 

buffer at input ports 

(b) 8-bit latch at output 2 

ports 


(c) 8-bit buffer for control : 
line 8 


5) 7400 

6 ) 7432 

7) 7490 

8) 7493 

9) 74154 

10) 7416? 

1 1 ) 740901 
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(d) 8-bit address latch/ 1 

output data buffer 

Two HAND gates used in the 1 

decoder logic and two inverting 
gates used for strobing output 
latches 

Two OR gates used in the 1 

decoder logic 

Decade counter. Two of them have 2 
been used in cascade to get the 
interrupt clock from 8156 timer 
4-bit binary counter. Two of them 12 

have been used to make an 8-bit 
counter in the input port 
I/O decoder 1 

4-bit BCD PEM.Two of them have 8 

been used for 8-bit BCD PRM In 
the output port 

Hex inverter. Two inverting 2 

gates have been used in deriving 
the ’clear’ pulse from the input 
port ’read’ line. Thus one chip 
serves three ports 
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12) 5MHz X-tal 

13) 1 nf cer disc 
cap 

14 ) 6.8K res. 

0.1 pf cer disc 
cap 

15) 1000 pf 
electrolytic 

cap 


Crystal for on the chip clock 
generator of 8085 
Used in generating the 50 psec 
clear pulses to the counters in 
input ports 
Decoupling capacitors 
(distributed on the card) 
Filtering capacitor at 
the V QC pin 


1 

6 

6 

25 

1 


There are two V lines coming from the. card with a 
cc 

curreht requirement of about 1A on each. The PC layout of the 
card is given in fig. PCI (appendix). 

8.2.2 Cardiovascular system card 

The PC layout of the card is given in fig. PC2 (appendix). 
This card contains the circuit simulating the cardiac conduction 
system and the circulation system. The contraction of chambers 
are indicated by LED s fixed on the card and coming out of the 
panel. The circuit (alongwith component values) for circulation 
system and conduction system have been given in fig. 6.2 and 
fig. 6.3 respectively. The current requirement of the card 
is about 100 mA. 
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8.2,3 PUS card 

The PC layout of the card is given in fig. PC 3 (appendix). 
There are six sensory channels forming the afferent division of 
PNS. Por observational convenience and flexibility the receptor 
and action potential blocks have been made separately and the 
external patch chord connection is suggested by the dotted lines. 

HI - R4 are the receptor cells or the generator potential 
blocks. The circuit has been shown in fig. 3.6 . They 
incorporate the features of adaptation and spatial and temporal 
summation. The component values (ref fig. 3.6) are given in 
table 8.2 , 

Table 8.2 

Component values in the receptor cells 


Cells No. of excitatory No. of inhibitory Components 
inputs inputs 


R1 1 0 3^-1 80E, C.,-10 \x£ 

R 2 -330K, Rj- 1.5 K 
R 4 -15K, C-10 iif 

Q - Be 147 

R2 1 0 R 1 -180K, C.,-10 jif 

R 2 ~680K, R^ - IK 
R 4 - 15K, C-10 \if 
Q - BC 147 




R1-180K, 0 1 — 10 [if 
R 2 -680]I, Rj-1 K 
R 4 -15K, 0-10 *if 

Q - BC 147 

IL,‘-180K, C^-lO^jf 

R 2 ’-680K, R 3 ’-1 K 

R - 100K, Q*~BC157 

Inverting gate - - 
74C901/6 

# 

R r 150K, 0^10 \if 
R2-680K, R 3 -1.8K 
r 4 -22k:, C- 10 \if 
Q - BC 147 

ly-isoK, c 1 *-ionf 

R 2 »-680K, R 3 »-1.8K 

R- 100K, Q’-BC 157 

Inverting gate- 
740901 /6 
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SI - S6 are action potential blocks with the circuit 
shown in fig* 3, 9. This circuit has also been provided on the 
card with the provisions of parameter variation* The panel shows 
the internal blocks of this circuit. The resting threshold, 
threshold decay rate, pulse width and absolute refractory period 
can be controlled using variable resistances provided on the 
panel. We refer this block as S. There is a simulation of a 
spinal reflex incorporating a synapse and a motoneuron. The 
motoneuron also has been simulated by the circuit in fig. 3.6. 

The component values for SI -S6, S and MN the motoneuron in 
spinal reflex are given in table S.3. 

Table 8.3 

Component values in S1_S6,S,MN 


Component 

R 

(K) 

R1 

(K) 

R2 

(K> 

R3 

(K) 

R4 

(K) 

EB 

OO 

01 

(pf) 

02 03 Q 

(pf)( f) 

Inverting 

gates 

SI - S4 

100 

51 

51 

15 

00 

100 

10K 

IdK .0.47 BC157 

740901 

S5,S6 

100 

51 

51 

6.8 

00 

100 

10K 

i 

10£ 0.47 BC157 

74C901 

S 

100 

6.2- 

220 

-6.2- 

220 

-15- 

82 

18- 

oo 

■ 100 

10K 

10K 0.22 BC157 74C901 

MN 

100 

51 

51 

15 

00 

100 

10K 

10K 0.47 BC157 740901 



86 


Measurement results for SI, S5 and S are presented in 
table S.4 and pulse repitition frequency YS input voltage plot 
for SI is given in fig. 8.2 . 


fable 8.4 

Measurements on SI, S5 and S 

(A) 11,22 


Block ) 

Si 

S5 

S 


11 (ms) 

0.59 

0.50 

0.1 - 

2.3 

T2(ms) 

0.45 

0.45 

0.5 - 

2.5 

lime period f 

and prf f for Si 

and S5 


Input voltage Y^ 

(v) 

3? (Si ) 

(ms) 

f (SI) 

(Hz) 

2(S5) 

(as) 

f(S5) 

(Hz) 

0.5 

27.6 

36.2 

9.4 

106.3 

1.0 

15.0 

66.6 

6.5 

153,8 

1.5 

11.4 

87.7 

5.0 

200 

2.0 

8.3 

120.5 

3.9 

256 

2.5 

6.3 

153.7 

3.1 

322.6 

3.0 

4.6 

217.3 

2.3 

434.8 

3.5 

3.4 

294.1 

1.9 

526.3 

4.0 

2.2 

454.5 

1.3 

769.2 

4.4 

1.73 

578.0 

1.1 

909.0 
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Ihe circuit for the synapse in the spinal reflex is the 
same as that for receptor in fig. 3.6 -with only one excitatory 
input and no adaptation. The component values are given "below 

HI - 22K, R2 - oo, R3 - 6.8K, R4 - 220K 

Cl - 0, C-10 \if Q - SCI 47 

The outputs of the six action potential blocks may be 
used for providing input to the six input channels of CNS card. 
The four output channels coming from CHS card are buffered by 
the hex buffer T4C902 to give a 0-5 V swing. Two outputs may be 
used for controlling cardiac activity by making the suggested 
patch card connections. The other two may be connected to muscle 
cell blocks (circuit shown in fig. 7.3). 

A reservoir of oxygen and other vital chemicals is 
simulated as shown in fig. 6.4. The output of R5 gives the 
level of reserves. The variable resistance simulating the 
physical exertion is available on the panel. R6, as shown in 
fig. 8.3 is a buffer for sensing the pressure in the circulation 
system. 

8.3 Operations and Experiments 

The system offers the possibility of carrying out a 
variety of experiments and investigations. For the convenience 
in carrying out these experiments, a description of the various 
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labeled 
GP 1 - 


GP 

AP , API 

AP 7 

AP 8 
AP 9 
AP 10 
AP 11 
AP 12 
AP 13 
AP 14 

AP 15 

MF 1 
ME 2 


pins on the panel is being given below 

<2?6 - Graded generator potentials from receptors 

R1 - R6. They may be used as inputs to action 
potential blocks. 

Graded voltage signal variable by a potentiometer. 
This may be used as input for studying the 
adjustable parameter action potential block (S). 

-6 Pulse outputs from the action potential blocks 
S, S1-S6. They may be used as inputs to CIS. 

Pulse output on the spinal reflex motoneuron. 

It may be applied to the skelectal muscle blocks. 
Pulse output on the output channel Ml of CIS, 

Pulse output on the output channel M2 of CIS. 

Pulse output on the output channel Cl of CIS. 

Pulse output on the output channel C2 of CIS. 

SA node output 

AV node output after the common His bundle. 

Action potential developed in the muscles of 
the right ventricle. 

Action potential developed in the muscles of 
the left ventricle. 

Force developed in skeletal muscle SMI. 

Force developed in skeletal muscle SM2. 
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MP 3 Porce developed in the muscles of right 

atrium (RAM). 

MP 4 Porce developed in the muscles of right 

ventricle (KVM). 

MP 5 Porce developed in the muscles of left 

atrium (1AM). 

MP 6 Porce developed in the muscles of left - 

ventricle (LVM) . 

HP 1 Blood pressure in the right airrium. 

BP 2 Blood pressure in the right ventricle. 

BP 3 Blood pressure in the lunges. 

BP 4 Blood pressure in the left atrium. 

BP 5 Blood pressure in the left ventricle. 

BP 6 Blood pressure in the aorta. 

BP 7 - BP11 Blood pressures in the downstream systemic 
circulation. 

The cardiovascular system can he studied separately hy 
leaving the patch chord connections to the excitatory and 
inhibitory control inputs of the SA node open. How the system 
works under the control of its own automatism. A delay in the 
conduction of excitation from SA node to ventricles may he 
observed. Por simulating an SA block, the conduction path from 
the SA node to the AV node has to be left open. How both the 
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atria and the ventricles -work at the AY node frequency. As 
the delays from the atria to the ventricles are nearly the 
same, there is no delay between the contractions of atria and 
ventricles, for simulating an AY block the conduction path , 
between the AY node and the ventricles are to be removed. A 
conduction disorder corresponding to the block in the common 
His bundle branch (section Q. in fig. 5.3) may be simulated by 
connecting the inputs to the ventricles together. In this case, 
the two ventricles contract in synchronism, but unrelated to 
the atria. His bundle branch block may be simulated by breaking 
the conduction path from AY node to one or both the ventricles. 
The ventricle not receiving the pulses through the AY node runs 
on its own at a relatively low frequency. Thus the conduction 
system simulation may be used for showing the various pacemaker 
levels and complete blocks of different kinds. However as the 
impairment in the conduction path is being achieved by breaking 
the path, partial blocks can not be simulated. The effect of 
external pacemakers may be studied by introducing pulses 
(TT1 levels) from an external pulse generator in the conduction 
path. 

Effects of controls from CHS on the SA node frequency 
may be studied by introducing pulses at the excitatory and 
inhibitory inputs of the SA node. These may he derived either 
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from two outputs of CHS or from external sources. 

Results of measurements on the conduction system is 
given below 
SA node 

(i) Frequency with both control inputs open 72 bpm. 

(ii) Freque cy with excitatory input of 500 Hz 

(pulsewidth-0.5 ms, 0r5Y swing) and inhibitory 
input open - 90 bpm. 

(iii) Frequency with inhibitory input of 500 Hz pulse 
width - 0.5 ms, 0-5V swing) and excitatory input 
open -66 bpm. 

AV node 

Free running frequency - 53 bpm. 
left ventricle oscillator 

Free running frequency - 41 
Right ventricle oscillator 

Free running frequency - 36 

Upper limit of frequency of excitations in the conduction 
system - 110 bpm. 

Belay from the SA node to the ventricles 100 ms 

However, the oscillator frequencies may be varied over 
a wide range by changing the settings of preset -potentiometers 
providing the bias to them. 


Pulse widths 


100 ms 
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A failure in the mechanical contractions of any of the 
ventricles may be simulated by breaking the connection between 
that ventricular oscillator and the corresponding muscle. 

Effects of transfusion an bleeding on the circulation 
system may be studied by injecting in or leaking out charge 
from the circulation system, for this purpose two free 
potentiometers provided on the left side of the panel may be 
used. Effects of physical exertion on oxygen and vital chemical 
reserves may be simulated by the potentiometer on the top as 
indicated on the panel. 

The inputs to the sensory neurons may be provided by the 
potentiometers on the left side. There is the provision of 
providing external sL gnals through the jack inputs as shown in 
fig. 8.4. These signals may be derived from some appropriate 
transducers. However, they should be in 0-5^ range. The 
potentiometer may be used to control the magnitude of the signal. 
In case, one is not interested in the generator potential block, 
the signals from the potentiometer may be directly given to the 
action potential blocks. The efferent outputs from the CHS 
may be given to the skeletal muscles shown on the panel and the 
SA node controls. Alternatively, they may be used for driving 
external actuators with the help of appropriate interfacing. 
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OHS card is available on the panel through input and 
output channels only. However, there is complete flexibility 
as regards the choice of connections. A possible set of 
connections is suggested on the panel. Various kinds of 
conditioned reflexes, learning processes and autonomic control 
activities may be simulated by appropriate software stored in 
the ROM on the CNS card. 

In writing the program, one has to consider the 
following factors - 

(a) Cycle time (T) - The time interval at which the processor 
is interrupted. Execution time for the program must be 
within cycle time, 

(b) Input frequency (f^) - Each input port is read at the 
start of a new cycle time, The count Hh accumulated in 
the counter is given as 

N. = f. T (8.1) 

To avoid a possible overflow, must be less than 

256. 

(c) Output frequency (f ) - The data to each output port is 
refreshed once in each cycle time. The 8-bit PEM in the 
output port is a BOD one and hence appropriate conversion 
should be done if the results after the processing are in 
the binary form. The timer output from 8156 serves as the 
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clock input to the PRMs. This is divided by 100 to provide 
the interrupt on RST 7.5 line of 8085. Thus the frequency of 
the PKM clock is given as 

f PEM = ( 8 * 2 ) 

If the data latched on the PBM input he Ho (BOD) then 
the output frequency is given as 



(8.5) 


8,4 A Program for Testing the Hardware of CHS Card 

lasting of a program for testing the hardware of CHS 
card is given in the appendix (TST). It follows the program 
structure suggested in section 7.1.5. 


Por a cycle time of 100 msec, the timer output frequency 
should he 1 KHz which requires a divi sion by 2500 of the system 
clock (2.5 MHz). A flow chart description of the program is 
given in fig. 8.5. Trap is kept low. After application of a 
RESET pulse, a frequency of 1 KHz from the timer indicates that the 
organisation of 8085, 8156 and 2716 is all right. A proper 
matching of input port and output port frequencies means that 
the entire hardware organisation of the card is in order. 
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Pig. 8.5 Plow diagram for testing program (1ST) 
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8*5 A Program Simulating Simple Conditioned Reflexes and 
Autonomic Control of SA node 

In modelling a simple conditioned reflex 

(fully deterministic behaviour), we take the case of a response 
R to a direct stimulus and an associative stimulus S a * C is 
a coefficient of association between S d and S Q and is refreshed 
in each cycle according to the relation 


C 


i+1 


°i + A 


0 . 


1 


(8*4) 


where Z> 0^ is a quantity depending on the stimuli as shown in 
fig. 8.6(a). The response R is determined by S a , S & and C as 
in the manner shown in fi- g» 8.6(b). 

In present studies, Si, S2 and Cl have been taken as 

« o d rn _ i = T , The conditioned reflex is 

S d’ a ana it * sd sa s 

modelled by the simple algorithm presented in fig. 8.6(c). 

Another conditioned reflex is modelled with S4 as the 
direct stimulus, S2 and S3 as indirect stimuli and M2 as the 
response. The manner in which C is transformed, the response 
and the algorithm for simulating this reflex are shown in 

fig. 8.7. 

Por the autonomic control, oxygen level is sensed through 
r 5 and S 5 * Pressure is sensed through R g and S g . C 1 and C 2 are 



I 



(c) Algorithm 


Fig. 8.6 Conditioned reflex involving S1,S2 and Ml 
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(a) Dependence of C* 0 n S2, S3, S4 
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SIM 





Fig. 8.9 Flow diagram of fee program SIM 
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A RESET pulse (by pressing RESTARE on the panel) dears 
all the past experience and the system starts learning afresh. 

8.6 Suggestions for Experiments and Further Work: 

The system developed offers flexibility for simulating 
a variety of structures and behaviours. By coupling transducer 
signals at the input to the nervous system and replacing muscle 
simulation circuits by appropriate actuators, a realistic 
simulation for many activities of living organisms may be 
developed. 

In the simulation of simple conditioned reflexes 
described in the previous section, the duration of ’learning* 
and 'forgetting* may be increased by using multi byte word to 
represent C. As only a small fraction of the cycle time is 
used for the execution of program simulating conditioned 
reflexes, one may write software for simulating much more 
involved learning processes with multiple loops of conditioned 
and unconditioned reflexes-making use of full cycle time. 

It is possible to use the input potentiometers (or 
whatever signal sources are being used) to provide certain 
parameters to the program. One possible way is making -.use of 
the TRAP. When TRAP is pressed, the program records the 
various input port readings and interprets them to give these 
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parameters* In an alternative approach, the input readings 
in the cycles just after the application of BESET pulse may he 
used to define these parameters. After defining these 
parameters, the inputs values are accepted as stimuli. By using 
these techniques, certain degree of flexibility in the program 
is achieved as one can vary the parameters like learning speeds, 
stimulus thresholds etc without changing the program stored in 
the EPROM. 

The work can he further extended in several directions. 

One may develop more elaborate models for simulating the neurons. 
In case miniaturisation of the circuit be possible, these 
models may be used for networks simulating reflexes and possibly 
many involved processes in CES. The simulation of the 
cardiovascular system suggests many possibilities for developing 
much more elaborate simulations. In place of treating pacemakers 
as localised oscillators, a simulation for distributed activity 
of each pacemaker may be simulated. Then the partial blocks 
in the conduction system and fibrillation etc also can be 
simulated. The circulation system simulation may be made more 
realistic by using a large no. of muscle elements in the 
musculature of heart chambers and thus simulating a propagation 
of contraction wave and making the circulation system much more 
distributed. Furthermore, linear resistances in the circulation 
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path may be replaced by nonlinear resistances simulating the 
constriction and dilation in blood vessels. 

By choosing a faster processor in the CITS simulation 
the no. of input and output channels can be increased enabling 
the simulation of more realistic processes. Sy introducing 
a random quantity statistically defined in a certain way m 
addition to C in relation (8.4) for transforming 0, random 
processes in the reflexes may be simulated. 
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APPENDIX' 

A. Linear Applications of CMOS Inverter 

The input~output characteristic 
of a CMOS inverter gate with = 5V y 
is shown in fig. A. 1 . It has a very 
■ high.. gain, .at- -VI •= ..-and- -a. high---injw^- 

impedance. It can be looked upoh as 
an Op Amp with the +ve input internally Fig. A. 1 
returned to and the -ve input being accessible from outside. 
This gate can be used as an Op <&np for many applications except 
of course as a voltage follower, because +ve input is 
unavailable. 

Rf 

Inverting Summer 

An inverting summer is shown 

in fig. A. 2 . From KCL at the input V 

.V 

v.-v V -V 

i x . o x A 

2 + ^7^ = 0 Kg. A- 2 






Ill 


For a single input 


V rt = V . - 
o t 


5t 

E i 


( T t - T t ) 


for two inputs ; 


V o “ V t 


H. R - 


R 


f ( \ - T t )- (V T t J 


Comparator 

Making use of one unity gain inverting circuit and 
another gate for summing, a comparator may be realised. 



= 


v( + v 2 W 27 t _ v +10. 

—5 — = 2 - t 2 


or V x - V t = 


V T i 


Hence for V 2 >V, ; V x -V t > 0 ^>Y o = »» 
for V 2 >V, ; Y x -V t > 0 =>Y 0 = High 

the + input and t’ 2 is the - input 


Thus V 1 is 
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For increasing the gain and reducing the error range, 
two more gate may he cascaded. The circuit is shown in fig.il. 4 . 





» /v v — J I 


y. 


Fig. A . 4 

Input V+ must be between 2V^_ and so that the inverter 
is in its linear range 


Error in the comparator 


R+ R 



Error in the comparator may arise from two sources: 
tolerance in the resistances and mismatch in the threshold 
voltage for the gates. Referring to fig. A*5. 


± 


Y 


tl 


'tl 


1 


R. 


R^ + Zi R^ 


A R. AR. 

< 2+ -57 > - < 1+ -17 > v i 

v x = (R 2 V / +(R 2 R 2 )Y 2 )/(2R 2 +AR 2 ) 
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=£> V - V. 


t2 


(1 + — - — ) V t1 - V 2 
2R 1 2R 2 1:1 ^ 


+ vzi + vfi 


2R 2 T 1 ^ R 1 
T~ 


error 


2 2R 2 ^ R, 

- — > 7 . 


1 


2R, 


1 


2R. 


tl 


V 2 + V 1 


£ 

2R, 


^ ZiR 1 
2 R„ 


^R 


Max error = Y^ + 


41 

R 


Y 


2 £SR 


T 


R 


V 1 

T 


= ^ + 


■ ^R 

R 




T i> 


Thus error has got two terms, one depending on Y'^ & Y ^ 
while other on z^Y, and V^. 

Differential imp 

A differential amp can he made with two CMOS inverter 
gates, with one as inverter and other as inverting summer as 
shown in fig. A. 6, 


R 1 R f 




or 


1 14 


V 1 ' = 2V t ~ V i 

v(- V. 7«-V T - V. V 

-2 — i + j-l + 1 + „2 — £ = 0 

R R R R 

V_v V - 7 Y - T 
V o= V t- R f(— 1 + ~ 1 + 1 ) 


R 


R 


R 


R. 


R 


£ (V, - V + 7 


V - V. 

t 

R 


R. 


If R f = R x ana V « 27 t 


Then 


V = 


R 


( v 1 - t 2 ) 


B; PG Card Layouts 

The PC card layouts for the three cards in the system 
are given in fig. PC1-PC3 as the following 

PCI (a) : CRS card (component side) 

PCI (b) : CBS card (back) 

PC2 (a) : Cardiovascular system card (component side) 
PC2 (b) i Cardiovascular system card (back) 

PC3 (a) : PNS card (component side) 

PC 3 (b) : PNS card (back) 


C, Program Listings 


Listing of the two programs TST and SIM are enclosed. 






PC 1 Cb)> Central Nervous System Ci 










.1 oU) Cnraiovanoal SralMft (mm***™* 












.A 1 1 ’“**“# 




PC J(b) Peripheral FsriroHs Sys'.nm fbaok’i 




BCD1 


RS75 


7 read all. in-ports 


ftCWJ ASSEMBLER, VfcR 2.0 ERRORS a u PAGE 


C 3.0001 

C32CU1 

SEC 4 
D30C 
3E49 
J 3 Q 0 
2 1 F 0 0 8 
r 9 

3 EOF 
0 308 
3 UR 
30 
FB 
00 

C3140) 


47 

AF 

0E08 

37 
7 8 
17 
47 
/ A 

022501 

14 
82 
27 
0 0 

C21C01 

Z9 

210008 

OBOO 

II ' 

DB01 

77 

23 

0b02 

n 

yfoa 


TSTs ORG OOOO'H 
JMP SIR i 
ORG OuSCH 
vJWP Rs75 
ORG 0100R 

strt: mvi A.'0C4ft 

OUT OCH 
MV I A.49H 
OUT 00 H 
LX I H t ObFOH 
SPHL 

M VI A , 0 C FVi 
OUT 08H 
mv i h,im 
DB 30B 
WAIT: El 
LOUP: NOP 

JpiP LOOP 


8CD2 : 


MOV B , A 

AFa a 
MV I C » 08H 
MOV 0,A 
MOV A,B 
«AL 

MOV B f A 
MOV AjrO 
JNC BC02 
INK D 
ADO D 
DA A 
OCK C 
OMZ BCDI 
' RET 

LX I H, 0800H 
IN OCR ' 

MOV M,A 

!»x U . 

Ik QIH 
MOV M,A 

XWX h' a 
U OSH 
MOV M, A 
XftX H 

IK U4H 
MOV M f A 


; LO*R k BYTE JF i 1 i G K 
,* HIGHER BYTE UP TIMER 
! SET SP 

; se r csr 

fSEX IN 4 R MASK 


; LOOP 

? BINARY to BCD 
MUA-BINARY 
;uHT: A-BCu 
;USES:B,C,0 




>'■ ACRj ASSEMBLER, VER 2.0 e. PROPS 


0 PAGE /. 


0 1 6 j 

0 i 6 A 

015.1 

01 so: 
01 6F 
01 7u 

0171 

0172 
0 1 7 6 
01 7o 
0 1 7 * 
0 1 7 A 
0 1 7b 
0 1 7 C 

017 u 
0 17 £ 

0 t 7F 

0 1 8 i 

018.1 
01 86 

018 7 


2 i 

/ E 

120708 
2 3 
7 E 
23 

BE 

F A 7 6 0 .1 
7 E 

320808 
23 
7 E 
23 
7 7 
2 1 
BE 

F 2 « i 0 t 
7 E 

320909 
f B 

’ Q 

w . 


PASO : 


PAS1 ; 


NO PROGRAM ERRORS 


If*X H 
MOV A | M 
ST A 0807H 
IRX H 
MOV A ,M 
I NX ' H 
CMP V 
Jm PASO 
MOV A.M 
ST A 0808R 
INX h 
MOV A , <* 

IN X B 
Mu V M , A 
I NX H 

cmp M 

Jr PAS1 

MOV A.M 

ST A 0809H 

El 

RET 

END 


U)UT2 = 1im2 


# u U 1 3=0 IN Up 3 , T * » ) 


?Q(IC4 = MaXIU5,I *b ) 


8080 

macro assembler. 

VER 2.0 

ERRORS to 

0 PAGE 3 



; 

SYMBOL 

TABLE 


* 01 

'' *’/ „ ' ' > 1 





A 

8 C 02 

H 

PASO 

SP 

0007 

0125 

mt 

00 m 

B 

C 

h 

PAS 1 

STRT 

0000 

8881 

0183 

0100 

' BCD 

k> 

TST 

0118 

Sill 

0006 

0000 ♦. 


•■mi 


asm 


H373 

WAIT 


8i,df •) 4 C R ;J ASSEMBLER, VaR 2.0 ePkORS a 0 PAGE 1 



:8 57 

•9 7 8 

A 17 
S 47 
C 7 A 
0 023 101 

9 H 


0133 00 

0134 C 2 2801 

013/ C*J 
0138 210908 

013d ?K 
013C C 92 401 

01 3F QiOA 

0141 2d 

0142 7E 

0 143 C 02 4(4 


JNZ BC01 
jREf . 

•XI H, 0809H 
IUV A/M 
ALB BCD 


0 U 0 y 
> u 0 0 

00 3C 

0 0 3 C 

0 1 0 0 

oi oo 

0103 

0 i 0 4 
0105 

01 Ob 
010A 

0 a o a 

01 oc 

OlOD 

0 i I o 
0112 
0114 
0 1 1 1> 
Oil 8 
0114 

one 

01 IE 
01 1 F 
0120 
0121 
0124 


Z 3 0 0 0 1 

C 3 3801 

2 J POO 8 

F0 

AF 

210008 
06E0 
7 7 

2 3 
05 

C20A01 

3 EC 4 
0 3 0 C 
3 £4 9 
D3Q0 
3FCF 
0308 
3d 1 B 
30 
FB 
00 

C32001 


SIM :OBG 0 0 0 0 H 
JMP STB 1 
OKG 003CH 
JMP PS 7 5 
„ OKG 0100H 

STRT: LX I H,0«F0H 
, SPHL 

IN IT ; XRA A 

LX I HfQROQH 
_ MVI o , OEOK 

IOITI ; MOV M, A 
I NX h 
DCH B 
JNZ IN1T1 
TIMM; MV I A . 0 C 4 rt 
OUT OCH 
MV I A r 49H 
OUT OOH 
; MVT A.OCFH 
OUT OoH 
MVI A-1BH 


! SF 1 *SP 

r iwi riAiBSrt miH 


; 1 IMFK LOWER oY f £ 

? 1 TMFK iiTGHgP BIT h 
7 SET CSK 


OE08 


WAIT: 

Db 3 OH 

? SF 1 1 M J’U 

El 


loop: 

NOP 


a, oop 

BCD: 

jmp 

LOOP 



IBIUAHV TO BCJ 
MNTA-SINAP’f 
* OUT: A-dCD 



BCDZ: 



;uS£S:b,c,d 

MOV 

B, k 



XRA 

A 


BCD! : 

MVI 

C , 0 8 B 


MOV 

0,4 



MOV 

At B 



7 WHITE OUTPUT PORTS 




A C«u ASSEMBLER 


ERRORS * 0 PAGE 


'JlbJ OB02 
U 1 6 b 86 
016b i f 
0167 77 

Olbti 26 
‘71 6y DBO 1 
0160 86 
016C IF 
716 D 77 

Oldfc 28 

Oi 6F dhoo 

oi n ee 

OH 2 IF 

0173 77 

0174 

on 4 1 A U 4 y H 

H? 0619 

0179 FAB 101 
one j A 0 5 0 8 
01 7 F D614 

0181 70 

0182 47 

0 181 C 6 3 2 

0185 120808 

0180 78 

0189 2F 
018a 3C 

0186 C632 

0 1 8 u 120908 

0190 16 IE 

0192 1 E7F 

0194 2690 

019b 2 El 0 

019? Ia 0008 

lilt |«m 

01 A| ^AfOl 


1ft 02H 
ADO M 
RAH 

MUV M , A 
OCX H 
l« oi h 
ADD M 
PAR 

MOV M , A 
OCX H 
IN 0OH 
ADD M 
PAR 

MOV M, A 

HEART: 

EDA U804H 
SUl 25 
DM HEAR! 
LDA OBuSil 
SU1 20 
HEARl : nop 

MOV B f A 
ADI 50 
STA 0808H 
MOV A,6 
CMA 
INK A 
ADI 50 
STA 0809H 

CNDRF: mV 1 0,30 
MVI £,127 
MV I H,90H 
MVI L # 1 Oh 
CR10 : LDA 0800H 
CMP 0 
JM CPU 
LDA 0 8 0 1 H 
CMP ... D 

LDA C oioAB 

MR A 
anz cri 7 

CRH : fl fllfl 
CR 12 ? LDA oUah 


cardiac C 7 . r 

COMP OXY r RCF 
DIFF=0XY-REF 

IF CUY>R£F C . 
Ulr Fa p j i E S S *• R F 
D I F 4 ’ = ? . Hi Fp 
SAVE 0 IFF 
XNh= 1 'J fi 0 + 1) 1 F F 


;exc*exco-diff 

CONDITIONAL REFLEX 
INPUT THRESHOLD 
ASSOC X AT 10tf» THRESHOLD 
; HIGH OUTPUT 
; LOW OUTPUT 




0 MACRO ASSEMBLER 


ERRORS s 0 PAGE 


»{?> FA3F02 
0 1 P to 3A0808 

o*p h ic 

,:'WFC C2Q1 02 

i, Ml t f If FF 

:m\ 420308 CR27 

iWj't 23 if 02 
MfJ 3A02Q8 CR22 

1 0,20 a y A 

f|03 f 23F02 

fin 060i 

|I| 3 021 102 

"i ||1 6 nF 

|21 / 320808 CR26 

02, U C33F02 

V Oglu 3AU108 CH21 

f f 2 Q a A 

022} FA3802 

, 3 A0208 

■022« f A3B02 

i A0 60 8 

02 2 C 0601 

| {jjj30 023402 

l ; |3f 320B08 CR28 

'« 023/ 68 

0238 F23F02 

’ *'0238 70 CP24 

02 3C C 34002 

§0708 gg|! 

83i? 51 

" ft 

40 PROGRAM ERRORS 


«■» * * a # 4S 3 3 

STA 08UBH 
JMP CR23 
LDft, 08023 
CMP 0 
JP CR23 
1»DA 080RH 
SUI 1 
JNC CF26 
XRA A 
STA oBOBil 
BMP CR23 
LDA 0801.} 
C*P 0 
JM CR24 
LDa 080 2 H 
CMP 0 
JM CR24 
LDA 080BH 
SUI 1 
JNC CR28 
XRA A 
STA 080BH 
CMp t 
JP CR23 
MOV A,L 


BOBO MACRO ASSEMBLER, VER 2.0 ERRORS » 


SYMBOL TABLE 


A 70519 


B£ - U % i - m -Pam - -S r m. 



